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Abstract 

We propose the existence of a hidden or dark sector besides the standard model (SM) of 
particle physics, whose members (both fermionic and bosonic) obey a local SU(2)h gauge sym¬ 
metry while behaving like a singlet under the SM gauge group. However, the fermiomic fields of 
the dark sector also possess another global U(1)h symmetry which remains unbroken. The local 
SU(2)h invariance of the dark sector is broken spontaneously when a scalar field in this sector 
acquires a vacuum expectation value (VEV) and thereby generating masses to the dark gauge 
bosons and dark fermions charged under the SU(2)h. The lightest fermion in this dark SU(2)h 
sector can be a potential dark matter candidate. We first examine the viability of the model and 
constrain the model parameter space by theoretical constraints such as vacuum stability and by 
the experimental constraints such as PLANCK limit on relic density, LHC data, limits on spin 
independent scattering cross-section from dark matter direct search experiments etc. We then 
investigate the gamma rays from the pair annihilation of the proposed dark matter candidate 
at the Galactic Centre region. We also extend our calculations of gamma rays flux for the case 
of dwarf galaxies and compare the signatures of gamma rays obtained from these astrophysical 
sites. 
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1 Introduction 


Experimental observations led by WMAP [T] and PLANCK [2] satellites reveal that only about 
4% of our Universe is made up of ordinary baryonic matter and about 26.5% of it is constituted 
by the unknown nonluminous matter or dark matter. Dark matter (DM), is supposed to have a 
very weak interaction with the visible sector of the Universe and strong evidence of its presence 
come from gravitational probes only. However, the particle nature of DM and reason for its high 
longevity (stability) are still unexplained. Direct searches of DM are also performed by various DM 
direct search experiments namely CDMS [3]-[6], CoGent [7j, XenonlOO [8], LUX [9] etc. Although 
no convincing detection of DM has yet been reported by these direct search experiments, until 
present, LUX experiment provides most stringent bounds on DM spin independent elastic scattering 
cross-sections usi of the dark matter-nucleon with respect to its mass. 

Dark matter particles can get trapped inside massive astrophysical bodies (due to their enor¬ 
mous gravity) like Galactic Centre (GC), the solar core. This may happen if the velocities of 
DM particles fall below their escape velocities inside those massive bodies. When accumulated 
in considerable amount, these trapped DM particles may undergo pair annihilation and produce 
fermion-antifermion pairs, y-rays etc. These y-rays or fermions (positrons, neutrinos, antiprotons 
etc), if found to be emitted in excess amount from these sites which cannot be explained by known 
astrophysical processes, may possibly be due to dark matter annihilation inside these sites. The 
detection of such excess y-rays, positrons, neutrinos etc. thus provide valuable indirect signatures 
of the particle nature of dark matter. Besides the GC, the dwarf galaxies may also be rich is dark 
matter. The dwarf galaxies are a class of faint and small satellite galaxies of our Milky Way galaxy. 
The huge amount of DM content within these dwarf spherical galaxies (dsphs) is inferred from 
their mass to luminosity ratio (^). The ^ ratio for these galaxies are found to be much higher 
than what is expected from the estimation of their visible mass. The dark matter rich dsphs can 
also emit excess y-rays due to the pair annihilation of dark matter. Nine such dwarf galaxies have 
recently been discovered in addition to the previously discovered 15 dwarf satellite galaxies of Milky 
way. Among these satellite galaxies, the y-ray flux obtained from the Reticulum 2 (Ret2 in brief) 
dwarf galaxy shows an excess of y-rays in the y-ray energy range of 2 to 10 GeV [10] while the null 
results from other satellite galaxies provide stringent upper limits on the annihilation cross-sections 
of DM particles for its various possible annihilation channels imiiii. 

In the last few years, the analyses of Fermi-LAT publicly available data |13j by several groups 
m-m have confirmed the existence of a low energy (few GeV range) y-ray excess which appears 
to be emerging from the regions close to the centre of our Milky way galaxy. The analysis of 
Fermi-LAT data m by Daylan et. al. [21] shows that the y-ray excess from the GC can be 
well explained by the annihilation dark matter scenario. They have also excluded all the known 
astrophysical processes which can act as the possible origin of this phenomena. In Ref. [21] it is 
shown that the observed y-ray spectrum from the GC can be well fitted by an annihilating dark 
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matter particle having mass in the range ~ 30 — 40 GeV which annihilates significantly into bb final 
state with an annihilation cross-section ~ (1.4 — 2.0) x 10“^®cm^/s with local dark matter 

density Pq = 0.3 GeV /cw?. In this work the authors have taken an angular region of 5^ around the 
centre of our galaxy as their region of interest (ROI) and used Navarro, Frenk and White (NFW) 
halo prohle with 7 = 1.26 for the computation of 7 -ray flux. However more recently, the authors 
Galore, Gholis and Weniger (CCW) of Ref. [22] have claimed to perform a detailed analysis of 
Fermi-LAT data along with all the possible systematic uncertainties using 60 galactic diffusion 
excess (GDE) models. The results obtained from the analysis of CCW provides a best ht for DM 
annihilation into bb final state having mass 49^54 GeV with (o'v)^j; = 1 . 761 q 27 x 10“^® cm^s“^. 
Moreover, the analysis of CCW assumes a NFW profile with 7 = 1.2 and a different region of 
interest (ROI) with galactic latitude |Z| < 20° and longitude \b\ < 20° masking out the inner 
region corresponding to \b\ < 2°. Different particle physics models for dark matter that are simple 
extensions of Standard Model (SM) are proposed to account for this 1-3 GeV excess in 7 -rays from 
GC [23| - |50| . It is to be noted that apart from dark matter, non-DM sources such as millisecond 
pulsars may provide a feasible explanation to the excess of 7 -ray observed at GC |5T|. Study of 
unresolved point sources near GC by Lee et. al |52| suggests that point sources also contribute 
significantly to the gamma ray excess. However, in this work, we will consider DM as the origin 
of the observed excess in GC gamma ray to explore the phenomenology of dark matter. Also the 
study of 7 -rays from previously known 15 different dwarf galaxies by Fermi-LAT m and eight 
newly discovered dwarf galaxies by Fermi-LAT with DES collaboration m give bound on DM 
mass and corresponding (av) for different annihilation channels. 

Since the Standard Model (SM) of particle physics can not possibly provide for the DM candi¬ 
date, an extension of the SM is called for. Different particle physics models for DM such as singlet 
scalar, fermion, vector where simple extension of the SM with a scalar, fermion or vector have 
been studied extensively in literature |53) - m- Inert doublet model (IDM) [72]-[ 86 ] also provides a 
viable DM candidate where an extra Higgs doublet is considered along with the SM sector. Various 
extensions of IDM with an additional scalar are also presented in Refs. m [ 88 ] . Several other 
models involving two Higgs doublet model (THDM) accompanied by a singlet DM (scalar, fermion 
or vector by choice) are also pursued in Refs. [8^-[92]. All these models in Refs. [53]-[H2] are based 
on a common approach where DM particle is stabilised by assuming a discrete symmetry {Z 2 or 
Z 3 ) and thus direct interactions (vertex with odd number of DM particle such as decay term) with 
the SM fermions and gauge bosons are prohibited. Hence, DM can interact with the visible sector 
only through the exchange of Higgs or scalar bosons appearing in THDM. Also there are models 
with multi component DM that presume discrete symmetry {Z 2 x Zg) [93] in order to stabilise the 
DM candidates. Different dark matter models with continuous symmetry such as U(l) or SU(2) 
gauge symmetry are also explored in literature [43], [Mj- |in4| . 

In this work, we consider a “hidden sector” framework of dark matter without pretending any 
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such discrete symmetry associated with it. We propose the existence of a hidden sector which has 
SU(2)h gauge structure. Dark fermions in this hidden sector are charged under this SU(2)h gauge 
group while all the SM particles behave like a singlet. Hence, the SM sector is decoupled from the 
dark sector and could interact only through the exchange of scalar bosons that exist in both the 
sector. Gauge bosons charged under SU(2)h are heavy and decay into dark fermions. Thus, the 
lightest one among dark fermions is stable and can be treated as a viable DM candidate. We show 
that the DM candidate in the present model that satisfy the limits from vacuum stability, LHC 
constraints, relic density, direct detection experiments can duely explain the Galactic centre 7 -ray 
excess and also is in agreement with the limits on DM annihilation cross-section obtained from the 
study of dwarf galaxies. 

The paper is organised as follows : In Sect. we present the hidden sector SU(2)h model. 
Constraints and bound on the model parameter space from vacuum stability, LHC results on SM 
Higgs, relic density etc. are described in Sect. In Sect. we show the results obtained for the 
available model parameter space and study of indirect searches of 7 -ray is performed. Finally in 
Sect. we summarise the work with concluding remarks. 

2 The Model 

We consider the existence of a “dark sector” that governs the particle candidate of dark matter. 
Just as the “visible sector” related to the known fundamental particles successfully explained by 
the Standard Model, we propose the existence of a hidden “dark sector” that relates the dark 
matter particles. We also presume that the Lagrangian of this hidden sector remains invariant 
under the transformations of a local SU(2)h as well as a global U(1)h gauge symmetries. Therefore 
we consider two fermion generations (* = 1 ) 2 ) where each generation consists of two fermions. 
Consequently, in the dark sector we have altogether four fermions namely fi (i = 1,4). The left 
handed component of each fermion (fn) transforms like a part of a doublet under SU(2)h while 
its right handed part behaves like a singlet under the same gauge group. Thus, the left handed 
components of /i, /2 and /s, /4 form two separate SU(2)h doublet^ However, both the left handed 
as well as the right handed fermionic components are charged under the postulated global U(1)h 
symmetry. The interactions between the dark sector fermions and the SM particles are possible by 
the presence of an SU(2)h scalar doublet through the gauge invariant interaction term 
which introduces a finite mixing between the SM Higgs boson and the neutral component of the 
hidden sector scalar doublet <h. This scalar doublet does not have any global U(1)h charge. As 
a result, the global U(1)h symmetry does not break spontaneously. However, being an SU(2)h 
doublet $ breaks the local SU(2)h symmetry spontaneously when its neutral component acquires 
vacuum expectation value (VEV) Vg. Besides the local SU(2)h gauge symmetry, the scalar doublet 

order to cancel the Witten anomaly | 1U6| we need at least two (even numbers) of left handed fermionic SU(2 )h 
doublets in our model. 
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which is in the fundamental representation of SU(2 )h gauge group, also possesses a custodial 
SO(3) symmetry. As a result of this residual SO(3) symmetry, three dark gauge bosons (z = 1 
to 3) which get mass due to the spontaneous breaking of the local SU(2 )h symmetry, become 
degenerate in mass. Non abelian nature of the SU(2 )h forbids the mixing between SM gauge 
bosons with dark gauge bosons (i = 1 to 3) [lOH 1105] . The scalar doublets H, <1> and the 
fermionic doublets can be written a^ 

’ ""“-“(/J,’ ’‘"‘■"I/J/ 

Therefore, the most general Lagrangian of the present proposed model contains the following gauge 
invariant terms 

£ D H - 

- As + Y, X^SV>'X^^) + E R) 

i=l,2 i=l,A 

-y'iXiL^/iR-2/2XiL^/2R>-y3X2L^/3R-l/4A2L^/4R + he , (2) 

with 

D, = + E + ’ 

a=l,3 

D', = + ( 3 ) 

a=l,3 

are the covariant derivatives of the SU(2)l x U(1)y doublet H and the SU(2)h doublets <h, 
respectively while <h = ia 2 ^* with (T 2 is the Pauli spin matrix. Moreover, g, g' and g^ are the 
respective gauge couplings corresponding to the gauge groups SU(2)l, U(1)y and SU(2)h. In the 
above equation (Eq. F'^^ is the field strength tensor for the gauge fields A'-^ (i = 1 to 3) of 
the SU(2)h gauge group while H is the usual SM Higgs doublet. The global U(1)h invariance of 
the dark sector Lagrangian forbids the presence of any Majorana type mass terms of the fermionic 
helds {fi, i = 1, 4) in Eq. We have assumed at the beginning that the dark sector fermions 
are charged under a global U(1)h symmetry. Therefore invariance of the dark sector Lagrangian 
(Eq. § under this U(1)h symmetry requires an equal and opposite U(1)h charges between each 
fermion and its antiparticle. Thus we can say that there is some conserved quantum number in 
the theory which can differentiate between a fermion and its antiparticle. In other words this can 
be stated the dark sector fermions in the present theory are Dirac type fermions. We have also 

^Although, in order to keep similarity with the expression of the Standard Model Higgs doublet H, we have 
introduced the notation of three scalar fields, in the expression of "I?, as G^, (f)^ and G®, however the symbols + and 
0 appearing in the superscript of dark sector scalar fields do not represent the electric charge of the corresponding 
scalar field as electric charge itself is not defied in the dark sector which is invariant only under SU(2 )h. 
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assumed that the dark sector fermions {fi, i = 1, 4) are in “mass basis” or “physical basis” so that 
the Lagrangian (Eq. does not contain any mixing term between these fermionic state^ The 
dark sector fermions can interact among themselves by exchanging dark gauge bosons and due 
to the presence of these interaction modes all the heavier fermions such as fi{i = 2 to 4) decay 
into the lightest one (/i). Consequently, the lightest fermion fi is stable and can be a viable dark 
matter candidate. Like the hidden sector gauge fields A'-^, the dark matter candidate /i also gets 
mass when the postulated SU(2 )h symmetry of the hidden sector breaks spontaneously by the VEV 
of 4>. Thus, the expression of mass of the fermionic dark matter candidate can easily be obtained 
using Eq. which is 


"i/i 


y'lVs 

^/2 


(4) 


We have already mentioned before, that due to the presence of the gauge invariant term 
the neutral components of both the scalar doublets, namely and (p^, possess mass 
mixing between themselves. The mass squared mixing matrix between these two real scalar fields 
are given by, 




scalar 


( 2Aiu^ 

X 3 VVS 

\ XzVVs 

2 X 2 v‘^ 


(5) 


After diagonalising the mass squared matrix we obtain two physical eigenstates hi and /12 

which are related to the old basis sates and (j)^ by an orthogonal transformation matrix 0{a) 
where a is the mixing angle between the resulting physical scalars. The relation between physical 
scalars hi and /i 2 with the scalar fields h^ and are given as 


hi = cos a h^ — sin a 


/i 2 = sin a h^ + cos a 


The expressions of the mixing angle a and the masses of the physical real scalars hi and /i 2 are 
given by 


a 


mi 


m 2 


1 _i 

- tan 

2 



\J Aiu2 + X2v‘l + - ^ 2 ^ 2)2 + [X-^VV^ , 

\l Aiu^ + X2v‘l - \f (Aiu^ - X 2 VIY + {X^WsY . 


( 6 ) 

(7) 


° Alternatively, one may think that the fermions in dark sector may have mixing between themselves similar to 
the case of SM fermions in quark and lepton sectors. Following the CKM mechanism in the qnark sector of SM we 
can assume that the mass matrix of up-type fermion generations i.e. fi and /s is diagonal while the mixing takes 
place between down-type fermionic states (/ 2 , fA- Now, since we have considerd the up-type fermion /i to be the 
lightest of all fermions in dark sector, thus in the present framework, the study of fermion mixing is redundant. 
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We assume the physical scalar hi is the SM-like Higgs boson which has been observed by the ATLAS 
and the CMS detector PUTHTUS] . Therefore we have adopted the mass (mi) of hi and VEV v of 
h^ to be 125.5 GeV and 246 GeV respectively. Thus, we have three unknown model parameters 
which control the interactions of the dark matter candidate fi in the early Universe, namely the 
mixing angle a, the mass (m 2 ) of the extra physical scalar boson /12 and more importantly, the 
mass mjj of the dark matter particle fi. In the rest of our work we have computed the allowed 
ranges of these model parameters using various theoretical, experimental as well as observational 
results. Throughout the work, for simplicity we take mass of fermionic DM candidate (/i) to be 
m. 

3 Constraints 

In this section we will discuss various constraints and bounds on model parameters that arise from 
both theoretical aspects and experimental observations. 

• Vacuum Stability - To ensure the stability of the vacuum, the scalar potential for the model 
must remain bounded from below. The quartic terms of the scalar potential is given as 

V 4 = + + ( 8 ) 

where H is the SM Higgs doublet and 4> is the hidden sector Higgs doublet. Gonditions for 
the vacuum stability in this framework is given as 

Ai > 0, A 2 > 0 , A 3 + 2\/ A 1 A 2 > 0 . (9) 

• LHC Phenomenology - In the present model of hidden sector (SU(2 )h) fermionic dark 
matter discussed earlier in Sect. an extra Higgs doublet is added to the SM. This dark 
SU(2 )h Higgs doublet provides an additional Higgs like scalar that mixes up with the SM 
Higgs. Large Hadron Gollider (LHG) performing the search of Higgs particle (ATLAS and 
GMS Gollaboration) have already discovered a Higgs like particle having mass about 125 GeV. 
The excess in 'y'y channel reported independents bv ATLAS jI07j and GMS m confirmed 
the existence of Higgs like bosons. In the case of Hidden sector SU(2 )h model, the mixing 
between SM Higgs with Dark Higgs results in two Higgs like scalars. In the the present 
scenario we take one of the scalar (hi) as the SM Higgs with mass mi = 125 GeV. We further 
assume that the signal strength of scalar hi also satisfies the limits on the same obtained for 
the newly discovered boson. Thus, hi in the present framework is identical with the SM like 
Higgs as reported by LHG Higgs search experiments (ATLAS and GMS). The signal strength 
of Higgs boson (h), decaying into a particular final state {xx, x is any SM particle), is defined 
as 

_ a{pp^h) Bvih^xx) 

aSM{pp^h)Br^^ih^xx) ' ^ ’ 
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where a{pp —?■ h) and Br(h —?■ xx) are the Higgs production cross-section and its branching 
ratio of any particular decay mode (x =quark, lepton or gauge boson), obtained from LHC 
experiments. The corresponding quantities computed using Standard Model of electroweak 
interaction are denoted by a^^{pp —)• h) and Br®'^(/i —)• xx) respectively. For the present 
model, the signal strength of the SM-like scalar hi is then defined as. 


= 


a{pp —)■ hi) Br(/ii —5- xx) 
{pp h) Br®^(/i xx) 


SM 


( 11 ) 


where the quantities are in the numerator of Eq. [m are the production cross-section and 
branching ratio of SM-like Higgs boson hi which are computed using the present formalism. 
Now due to the mixing of scalar bosons, the coupling of SM-like Higgs boson to the SM 
fermions and gauge bosons are modified with respect to SM Higgs boson [h) by the cosine 
of mixing angle a whereas, the couplings of non-SM scalar boson /12 to SM particles are 
multiplied by a factor sin a. Hence the ratio = cos^ a and from the similar argument 

one can yield = sin^a. The SM branching ratio can be expressed as Br®^(/i —)• 

xx) = -—where T^^{h —)■ xx) is the decay width of SM Higgs boson h into any 
final state particles and T®'^ is the total SM Higgs decay width having mass mi = 125 GeV. 
Similarly one can derive the expression for branching ratio of hi into any specific decay channel 
in the present model Br(/ii —>■ xx) = where ri(/ii —)> xx) = cos^ Q;r®^(/i —)> xx) is 

the decay width of hi decaying into xx final state while Ti is the total decay width of hi in 


the present model. Hence, the signal strength of hi in Eq. 11 can be written in the form 

pSM 


Ri = cl 


( 12 ) 


where we have denoted cos a as Cq. It is to be noted that apart from the decay into SM 
particles the SM-like scalar hi can also have invisible decay mode into dark matter particles. 
Therefore the total decay width of hi, in the present model, can be written as 


T-' _2 t^SM I T-'in 


(13) 


In Eq. 


13 


pinv 

^ 1 


is the invisible decay width hi for the channel hi —)• fifi- Eor mi > 2m the 


expression of invisible decay width of hi is given by 

4m2 \ 3/2 


pinv _mim 2 


1 - 


mt 


(14) 


since coupling between hi and dark matter candidate is proportional to ^Sa- In the above. 
Vs is the VEV of SU(2 )h Higgs doublet and Sa = sin a. Similarly for the other scalar 
involved in our model, the signal strength R 2 is expressed as 

a{pp —)• /12) Br(/i 2 


R2 = 


XX} 


(jSM(^pp Br®’^(/i xx) 


(15) 












with a{pp —)■ /i 2 ) being the production cross-section of /12 and Br(/i 2 —)• xx) is decay branching 
ratio of /i 2 to any final state. However in this case, the Standard Model predictions a^^{pp —)• 
h) and Br^'^(/i —)• xx) are computed for the mass of SM Higgs boson rrih = m 2 . Using the 
similar approach we used to calculate Ri and replacing hi, cos a etc. by /i 2 ,sina the signal 
strength R 2 of ^2 can be expressed as 


^T^^{mh = m2) 

U2 = 


(16) 


where = m 2 ) is the total decay width of SM Higgs boson if it has mass = m 2 

while r 2 is the total decay width for the non-SM scalar boson /12 


r 2 = slT^^imh = m 2 ) + rr 


(17) 


The coupling between dark matter and h 2 depends on the factor — coSq. Hence, invisible 

Vs 

decay width of h 2 (T^"'") for m 2 > 2m is given as 


r inv 
2 


m 2 m? 2 ( \ 

Stt ml J 


(18) 


As stated earlier, we consider hi with mass mi = 125 GeV to be the Higgs like scalar and infer 
Ri > 0.8 1 109) and invisible decav branching ratio Br;^.j, < 0.2 do] where Br(^^ = T^^^/Ti is 
defined as the ratio of invisible decay width to the total decay width. 


• Dark matter relic density - The DM relic density as measured by PLANCK satellite 



Figure 1: Feynman diagrams for dark matter annihilation into fermions (quarks and leptons), gauge 
bosons and scalars contributing to DM annihilation cross-section. 
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experiment is given as [2] 


Odm/i = 0.1199±0.0027 . 


(19) 


In Eq. 19 his the Hubble parameter measured in the unit of 100 km s“^ Mpc“^. We calculate 
the relic density for the fermionic (SU(2 )h) dark matter candidate in the assumed dark sector 
in our model by solving the Boltzmann equation. Relic density of the DM candidate is 
obtained by solving the Boltzmann equation m 


dn 

dt 


+ 3Hn = — (av)— nlq) , 


( 20 ) 


where n is the number density of DM particle and Ueq is the same in equilibrium. In Eq. 20 


(cTv) is the thermal averaged annihilation cross-section of DM particle into SM sector and H 


is Hubble parameter. Solution to Eq. 20 gives the DM relic abundance of the form 

1.07 X lO^Xf 


f^DMh = 


^/^Mpi(f7v) 


( 21 ) 


where is the effective number of d.o.f (degrees of freedom), Mpi is the PLANCK mass 
(~ 1.22 X 10^® GeV) and Xf = m/Tj with Tj being the freeze out temperature of the DM 
species respectively. We compute the freeze out temperature 7/ by solving iteratively the 
following equation 




= In 



45M2 

(cjv) 


( 22 ) 


In order to obtain the freeze out temperature of DM and hence its relic density using Eqs. 21 


22 we need to calculate the thermal average of the product between total DM annihilation 
cross-section (u) and the relative velocity (v) of two annihilating DM particles. The expression 
for the thermally averaged DM annihilation cross-section into all possible final states is given 


as 


(o-v) = 


1 


Sm^Tf K^{m/Tf) 


poo 

/ ds a{s) {s — y/s Ki{y/s/Tf), 
J Am? 


(23) 


where the factors Ki, (i = 1, 2) are the modified Bessel functions and ^/s being the centre of 
mass energy. In the present formalism dark matter candidate /i can annihilate into the SM 
particles through s-channel processes mediated by the scalar bosons hi and / 12 . In the above 
Eq. 23, a{s) denotes the total annihilation cross-section of dark matter into all possible final 


states which are allowed by the Lagrangian given in Eq. Feynman diagrams for different 
annihilation channels of /i are shown in Fig. The expressions of av for different final state 
annihilation of dark matter into SM particles are derived from the Feynmann diagrams shown 
in Fig. The value of cjv obtained for DM annihilation into SM fermion and antifermion 
pairs (//) at the final state is of the form 
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avff = N, 




C 9 

vi 


Stt 


1 - 


4 mj 


3/2 


(s — 4 m^) 


+ 


2 (s — mf)(s — m^) + 2mim2rir2 


(s — mf)^ 4-mfrf (s — 

( 24 ) 


[{s — m\Y + m\V‘f\[{s — m\Y + m\V‘^\ ’ 

where m is the DM mass and m/ is the mass of specific fermion (/ = quark or lepton). 
In Eq. ^ V and Vg are the vacuum expectation values of SM Higgs doublet and dark Higgs 
doublet, Nc is the colour quantum number (3 for quarks and 1 for leptons). Moreover, Ei, 
r2 are the total decay widths of the scalar bosons hi, h2 and the expressions of Ei and E2 
are given in Eqs. 13 Idil We also calculate av for W~^W and ZZ channels which proceed 
through the s-channel exchange of scalar bosons hi, /12 (see Fig. [^. The expressions of 
(y^w^w- s-^cl cjv^z are furnished below 




and 


o^zz = 


v 1 Stts 


+ •- 


1_ 4r^ 


1/2 




1 + 


(5/2 - TO^) 

2m^ 


(s — 4 m^) 


2 (s — ■mf)(s — m|) + 2mim2rir2 


(s — + toIeI [(s — mf)^ + mjr^][(s — 771^)^ + rrijr|] 


a ^a. 

9 ? 167rs 


4m\ 

1 -^ 

s 


1/2 


/ 2 to 

V V 


1 + 


(s /2 - to|)- 
2777^ 


(s — 4 m^) 


(s — m\Y + 'm\T\ 
( 25 ) 

1 


1 


2 (s — m\){s — m^) + 2mim2T 1T2 
(s — + rn^2 [(■^ — 'm\Y + [(s — m\Y + mf E^] 


(s — mlY + m\T\ 

( 26 ) 


In the above, mw and mz denotes the respective masses of W and Z bosons. Annihilations 
of DM particles into scalar bosons hi and /i2 are also taken into account. The process of DM 
annihilation into scalars hi or /i2 is also scalar mediated, depends on scalar couplings between 
hi and /i2- The s-channel annihilation cross-section of fi annihilating into the pairs of hi and 
h2, calculated using fifi ^ hihi, 7 = 1,2 annihilation diagram, takes the following form 


^^hihi — 


m 


10775 v1 


1 - 


4 m\ 4 m\{m\ — 1 ) 


1/2 


(s — 4 m^) 


x 2 


(s — m\Y + mfrf 


^2 X 2 
^a2'211 


2SaCaAlllA21l((s - mf)(s - m|) -h 2TOiTO2rir2) 


(s-^ 




[(s — 777^)2 -I- 777 fr^][(s — 777^)2 -|- m\T‘^ 


( 27 ) 


and 


^^h2h2 — 


m 


IGtts v‘i 


, 4777n 

1 -^ 


4777^(7772 -1)\ 

S 2 J 


1/2 


(s — 4777^) 


«2 X 2 
^a'^122 


{s — + TTi^r^ 


^2 \2 
^a'^222 


2SQ,CaAi22A222((s “ 777f)(s - 777|) 2777i7772rir2 ) 


(s — 7772)^ -I- 7772 r 2 [(s — 777^)2 -I- 777 frf][(s — 777^)2 -|- m\T‘^ 


( 28 ) 


where, Xijk is the coupling for the vertex involving three scalar fields hihjh^. The expressions 
for the scalar couplings Am, A211, A122 and A222 are given in Appendix A. We calculate the 
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thermally averaged annihilation cross-section of the present DM candidate using Eqs. 23]|^ 


We then compute the freeze out temperature 7/ by solving Eq. 22 and finally obtain the relic 


density of /i at the present epoch from Eq. 21 


DM Direct Detection - Direct detection of DM particle is based on the scattering of the 
DM particle with the target nucleus of the detector material. Fermionic dark matter in the 
present model can undergo elastic scattering with the detector nucleus. This elastic scattering 
of the DM and the nucleus will transfer a recoil energy to the target nucleus which is then 
calibrated. From the non-observance of such elastic scattering events the direct detection 
experiments give the upper bound of elastic scattering cross-sections for different possible 
masses of dark matter. The scattering cross-section is expressed as cross-section per nucleon 
for enabling direct comparison of the results from different experiments. In the present model 
DM fermion of mass m can interact with the target nucleus through t-channel Higgs mediated 
processes through both hi and /i 2 - The spin-independent (SI) elastic scattering cross-section 
off the detector material normalised to per nucleon can be written as 


crsi = 


sin^ 2a m? o 
- 


dvr 


vt 


2 ^2 I '^P 


m 


1 


nir 


(29) 


where nir = 


m-\-mp 


^ is the reduced mass for the DM-nucleon system and \p | 66 ] is given in 


terms of the form factors /„, proton mass rup as 


_ rr^ 

Ap — 

V 




L q 


~ 1.3 X 10 


-3 


(30) 


Using Eqs. 29 30, we calculate the spin independent elastic scattering cross-section of the DM 
fermion off the nucleon and compare it with the experimental bounds from LUX [9]. 

Note that both DM annihilation cross-section and DM-nucleon scattering cross-section depend 
on an effective coupling ^e// = \'^SaCa\ (Eqs. 24 29). This effective coupling is a useful 


parameters to explain the dark matter phenomenology in the present framework. Further 
discussions on the effective coupling are given later in Sec. 


• DM Indirect Detection The existence of DM has now been well established from gravita¬ 
tional evidences in astrophysical scale. Indirect search of DM focuses on the non-gravitational 
search of DM candidate and explores the particle physics nature of DM. The astrophysical 
sites such as Galactic Centre (GC), dwarf galaxies etc. are of great interest since dark matter 
can be trapped and accumulate at GC due to the enormous gravity in the region of GC 
and the mass to luminosity ratio of dwarf galaxies indicate the presence of dark matter in 
large magnitude. These sites are suitable for indirect search of DM as DM particles trapped 
in these regions can undergo annihilation into various SM particles which can further pro¬ 
duce gamma rays, neutrinos etc. Thus any observed excess in the fluxes of 7 -ray, positron, 
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anti-proton from such sites can indicate DM annihilation processes in those sites if other 
astrophysical phenomena cannot explain the observed excess. Fermi-LAT ra searches for 
the excess emission of 7 -rays originating from GC and dwarf galaxies. Observation of the 
excess in le~ and p/p flux is performed by AMS-02 [113] experiment. In this Section we 
will study Fermi-LAT observed gamma ray flux results from the centre of Milky Way and 
surrounding dwarf spheroidal galaxies (dSphs). 

The expression for the differential 7 -ray flux obtained from a region of interest (ROI) subtends 
a solid angle dD centered at GC is given as 

dEdD “ Svrm^^ ^ ^ dE^ ' ^ ^ 


where {crv) j is the average thermal annihilation cross-section of DM particles annihilating 


into final state particle / and is the photon energy spectrum of DM annihilation into the 


dE. 


same. The factor J appearing in Eq. 31 is related to the quantity of dark matter present at 
the astrophysical site considered and is expressed in terms of dark matter density as 


J= / p‘^{r{s,6))ds . 


(32) 


/ los 


In Eq. 32 the line of sight (los) integral is performed over an angle 6 , is the angular aperture 
between the line connecting GC to the Earth and the direction of line of sight. In the above 
Eq. 32 r = \JrQ + — 2rQS cos 6 where r© = 8.5 kpc, is the distance to the Sun from GC. 

It is clear from the expression of Eq. 32 that value of J factor is dependent on the nature 


of the chosen p(r) factor i.e, DM halo density profile p{r). In the present work, we consider 
Navarro-Frenk-White (NFW) |114j halo profile. DM density distribution for the NFW halo 
profile is given as 


P{r) = Po 


jr/rs) 

(1 -h r/rs)^~'^ ’ 


(33) 


where = 20 kpc is the characteristic distance and po is normalised to local DM density i.e., 
Pq = 0.4 GeV cm“^ at a distance r© from GC. 

The analysis by Daylan et. al. [2T] of Fermi-LAT data suggests an excess in 7 -ray in the 
7 energy range of 2-3 GeV at GC. The same analysis demonstrates that this excess can be 
explained by the annihilation of 31-40 GeV DM into bb with (crv)^,^ = 1.4 —2.0 x 10“^® cm^s“^. 
In this work [ 2 T], inner galaxy gamma ray flux (5*^ from GC) is calibrated using NFW halo 
profile with 7 = 1.26 and local DM density p© = 0.3 GeV cm“^. In a recent work by Galore, 
Cholis and Weniger (CCW) [ 22 ] detailed analysis is performed for the GC 7 -rays along with 
the systematic uncertainties using 60 galactic diffusion excess (GDE) models. Results from 
CCW analysis provides a best fit for DM annihilation into bb having mass GeV with 
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{^^)bb ~ ^ 10“^® cm^s“^. However, CCW analysis of Galactic Centre excess (GCE) 

for gamma ray have also considered generalised NFW profile ( 7 = 1 . 2 , pq = 0.4 GeV cm“^) 
for a different region of interest (ROI) with galactic latitude |/| < 20*^ and longitude \b\ < 20*^ 
masking out inner \b\ < 2®. In another work P. Agrawal et. al. |115j reported that annihilation 
of heavier dark matter (upto 165 GeV for bh channel) can also explain the observed GCE in 
7 -ray when uncertainties in DM halo profile (NEW) and the J-factor are taken into account. 
However in the present work, we do not consider any such uncertainties in halo profiles or 
J vaules and use the canonical NEW halo profile used in CCW analysis. Using Eqs. 
we calculate the 7 -ray flux (in GeV cm“^ s sr“^) for the ROI described in CCW analysis for 
Eermi-LAT data. As mentioned earlier we consider for our calculations the NEW profile with 
7 = 1.2 and pq = 0.4 GeV cm“^. 


Apart from the GC region, dwarf galaxies of the Milky-Way galaxy are also of great signif¬ 
icance for indirect search of DM as these galaxies are supposed to be rich in dark matter. 
Recent analyses of 7 -ray fluxes from 15 Milky-Way dSphs reported by Eermi-LAT m provide 
a limit on DM mass and corresponding thermally averaged annihilation cross-section {crv)^ 
into different channels / (r and b). Eermi-LAT have used their 6 year data collected by Fermi 
Large area Telescope and performed an analysis for 15 dSphs using “pass -8 event level anal¬ 
ysis" (see m and references therein). In an another work [12] Eermi-LAT in collaboration 
along with Dark Energy Survey (DES) collaboration also provide similar bound on {crv)j 
where they include data for 8 new dSphs. For both the analysis presented in mm a canon¬ 
ical NFW halo profile (7 = 1) is considered, and the astrophysical J factors are measured 
over a solid angle AH = 2.4 x 10“^ sr with angular radius 0.5^. Independent searches carried 
out by Eermi-LAT m and DES-Fermi-LAT collaboration on 15 previously discovered and 8 
recently discovered different dSphs reported no significant excess in observed 7 -ray. Results 
from the DES dSphs [ 12 ] also predicts an upper bound to the observed 7 -ray energy flux with 
95% confidence limit (C.L.) for 8 newly found dSphs. Gamma ray flux for dwarf galaxies 
when integrated for an energy range extending over a region of solid angle AH is expressed 
as 


(dv) r^max 

- 9 —^ / dh/^ , 

StTITIdm j Smin ^-^7 


(34) 


where 4w- is the 7 -ray. The expression of flux presented in Eq. 34 is calculated for a single 


final state annihilation of DM. Hence, summation over different final channels is not needed. 


Form of J factor appearing in Eq. 34 is different from Eq. 32 and written as 


J = [ f p^{r{s,e))ds , (35) 

J AQ J\os 

calculated over a solid angle AH = 2.4 x 10“^ sr subtended by the ROI (0.5*^ angular radius) 
for NFW halo profile (7 = 1). The density distribution function for NFW profile with 7 = 1 
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is then 


p{r) = po 


r 


3 

s 


r{rs + r 


2 ’ 


(36) 


where is the NFW scale radius and po represents the characteristic density for the dSphs. 
In the case of Fermi-LAT analysis, J factors for different dSphs are adopted from Ref. m- 
We use values of J factor from [12] for computing gamma ray flux for 8 DES dSphs for the 
dark matter candidates in our model. However, it is to be noted that J factors for DES dSphs 
candidates are obtained assuming the point like dSphs instead of having spatial extension (as 
in the case of m) to avoid the uncertainties in halo profile arising from spatial extension. 
Calculation of gamma ray flux is also based on the assumption that the spectrum follows 
the conventional power law ~ As mentioned earlier, study of 15 dSphs by Fermi-LAT 
and 8 other dSphs by DES-Fermi-LAT collaboration found no significant excess in 7 -ray from 
these dwarf galaxies. However, a recent search on a newly discovered dwarf galaxy Reticulum 
2 (Ret2) in a work by Geringer-Sameth et. al [l0| has reported an excess in observed 7 -ray 
signal. In the present work, we calculate the 7 -ray flux for annihilation of hidden SU(2 )h 
fermionic dark matter into 7 -ray through different SM final states and explore whether the 
model can account for GCE in 7 -ray and also satisfies the bounds on gamma ray flux from 
dwarf satellite galaxies. 

As mentioned earlier, in the present model dark matter candidate (/i) is fermionic in nature 
and it interacts with the visible world (SM particles) through the exchange of two real scalar 
bosons hi and /i 2 . As a result the annihilation cross-sections of the DM dark candidate /i 
into the final states that composed of SM particles (mainly light quarks and leptons) are 
proportional to the square of relative velocity (v^) between the annihilating dark matter 
particles (p wave process). Now the averaged DM relative velocity is proportional to 
nsj-mi! with X = ^ is a dimensionless quantity and T being the temperature of the 
Universe. Hence, in our model, the thermally averaged annihilation cross-section used for 
computing DM relic density, at x ~ 20 — 30, is different from the annihilation cross-section 
(for X ~ 3 X 10® [116j - [TT7] I needed to calculate 7 -ray flux at the Galactic Gentre and dwarf 
galaxies. The latter quantity is velocity suppressed as the average DM relative velocity is 
~ 10“® when the annihilation of DM occurs at the GG. Among all the annihilation channels 
of /i, the annihilation mode /i/i —)• bb plays a significant role for the 7 -ray excess observed 
from GG and dwarfs satellite galaxies as it is the most dominant annihilation channel for the 
considered mass range of DM. In order to explain the GG gamma-excess by DM annihilation 
to bb, the annihilation cross-section should be ~ 1.76 ^q27 x 10“^® cm^/s [22|. Although in 
the present case, the thermally averaged annihilation cross-section for the bh annihilation is 
quite small, however the quantity (o'v)^,^ can be significantly enhanced using Breit-Wigner 
resonant enhancement mechanism mni-iiizi. Breit-Wigner enhancement occurs only when 
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the mass of the dark matter (m) is nearly equal to half of the mediator mass (in our case it is 
the mass of /i 2 ). Therefore, we have defined the mass of the hidden sector scalar boson (h 2 ) 
and the centre of mass energy ^/s in the following way 

771,2 = 4r?i,^(l — 5) and s = 4m^(l + z) , (37) 


where <5 < 0 represents the physical pole and 2 ; is the measure of excess centre of momentum 

for the bb annihilation channel, can now be 

written as 


energy scaled by 47 ?i^. In terms of z, Eq. 23 


4 X , _ , _ 

'\^"^)hh= / dza{z\iz^\^ zKx{2xyJ\^ z) 

Jo 


with the expression of a{z )bb is given by[| 

/ N gc y/z 
(^{z:)bb = ] 

and 


1 + z- 


2 x 3/2 


4r?T,2 1 + 2; 1^(2; + (5)2 _|_ 


9c = 


Nc { m cos a mb sin a 


16 TT V Vs 


(38) 


(39) 


(40) 


where 72 = r 2 being the total decay width of h 2 of mass m 2 - It is to be noted that the 
upper limit of the above integration should be 00 (see Eq. [23| ) , however the integrand becomes 
negligibly small when 2 : approaches to z^ff ~ max[4/x, 2|(5|] for d < 0 [37]. |117j . Using the 
above prescription, we calculate the thermally averaged annihilation cross-section (cTv)^j; of 
the dark matter candidate /i for GC and dwarf spheroidal galaxies. The actual values of 
{av)jyl, 72 and <5 for the two chosen bench mark points (BPl, BP2) are given in Table of 
Sec. 1^ We have found that for |(5| ~ 10“^ the annihilation cross-section (uv ),5 ~ 1.9 X 10-26 
cm^/s which can explain the excess of gamma ray flux in GCj^ 


4 Calculational procedures and Results 


In this section we present the computation of dark matter annihilation cross-sections as also the 
DM-nucleon elastic scattering cross-sections. They are required for the calculation of relic densities 
and the comparison of the latest DM scattering cross-section bound given by the LUX direct 
detection experiment. The invisible decay widths and signal strengths for the SM-like scalar is 
also calculated in order to constrain the model parameter space. The gamma ray flux are then 
computed within the framework of SU(2 )h fermionic dark matter for galactic centre as also for 
dwarf galaxies and the results are compared with the experimental analysis. 

^Since the Breit-Wigner enhancement occurs when m ~ 7712 / 2 , as a result only the term proportional to 
"^41 dominanatly contribute to the annihilation cross-section appearing in Eq. 


24 


® Similar results for Breit-Wigner enhancement of dark matter annihilatin cross-section have been reported in |37| . 
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4.1 Constraining the model parameter space 


The fermionic dark matter in the present model can annihilate through scalar mediated (hi and 
/ 12 ) s-channel processes. As mentioned in Sec. the model parameter space is first constrained 
by the vacunm stability conditions given in Eq. The signal strengths Ri and R 2 for the Higgs 


doublets hi (SM) and /12 (dark sector) are then compnted nsing Eqs. 19 23 With the chosen 
constraints on Ri {Ri > 0.8, Ref. |in9p the invisible decay branching ratio of SM-like Higgs Brj^^ 
is calcnlated and the parameter space is further constrained by LHC experiment limit of Brj^^ 
(Br;^.^, < 0.2 [TOj). The parameter space thus constrained is then used to compute the thermal 
averaged annihilation cross-section (uv) of the present fermionic dark matter candidate and the 


relic density is obtained by solving the Boltzmann equation (using Eqs. 20- 23). The annihilation 


cross-sections are computed with the calculated analytical formulae given in Eqs. 24||^ with two 
choices of VEV for 4* (dark Higgs donblet) namely Vg = 246 GeV and 500 GeV. In onr calculation we 
consider the mass mi of the SM-like Higgs boson hi to be 125 GeV. The calculation is performed for 
two values of the dark sector scalar /12 masses and they are m 2 = 100 GeV and 110 GeV. These relic 
densities are compared with the dark matter relic density given by PLANGK [2]. Thus PLANGK 
result further constrains the parameter space of our model. With this available parameter space we 
evaluate the dark matter-nucleon spin independent scattering cross-section (dsi) for the purpose 
of comparing our results with those given by the dark matter direct detection experiments such as 
LUX, XENONIOO etc. In this way we restrict onr model parameter space by different experimental 
results. 

In Eig. and Fig. we show the calculated values of usi with different DM mass in the 
present model where the conditions from vacnum stability, bound on SM Higgs signal strength and 
DM relic density results from PLANGK have been imposed. We first choose certain values of mi 
and m 2 and vary the couplings Aj, z = 1 to 3 (satisfying vacunm stability conditions given in Eq. 
for two different values of Vg which also constrain the mixing angle a throngh the Eq. Here we 
want to mention that we have varied Ai and A 2 in the range 0 to 0.2 with the valnes of both A’s are 
evenly spread within the considered range. Gonsequently the value of the parameter A 3 becomes 
fixed by the vacuum stability criteria given in Eq. which is also varied with equal interval in 
the range IA 3 I < 2t/Ai^. The model parameter space thns obtained is then fnrther constrained 
by imposing the conditions Ri > 0.8 and ^ from LHG results. Using this restricted 

model parameter space satisfying both vacuum stability and LHG bonnds, we therefore calculate 


the relic density of the dark matter candidate fi by solving the Boltzmann equation (Eq. 20) for 


different values of DM mass. Finally, we consider specific range of model parameter space which 
is in agreement with DM relic density reported by PLANGK experiment and for these parameter 


space we compute the spin-independent direct detection cross-section using Eqs. ^ In this 
way the viable model parameter space for the dark matter candidate fi is obtained. Fig. is 
for the case m 2 = 100 GeV while Fig. is for the case m 2 = 110 GeV. The upper limit on <731 
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(a) (b) 

Figure 2: The allowed range of m — cjsi parameter space obtained for m 2 = 100 GeV (left panel) 
and m 2 = 110 GeV (right panel) plotted using the bounds from vacuum stability, LHG constraints 
on SM Higgs and relic abundance of DM obtained from PLANGK [2]. Limits on DM-nucleon 
cross-section from LUX [9] is also plotted (blue line) for comparison. 


for different values of DM mass, obtained from LUX DM direct search experiment, are also shown 
in Fig. [^-b by the blue line for comparison. The red and green scattered regions as shown in 
Fig. [^-b correspond to two choices of Us=246 GeV and 500 GeV respectively. From Fig. it can 
be observed that only the region near the resonances of scalar bosons hi and /i 2 is in agreement 
with the upper limit on usi predicted by LUX. It is also seen from Fig. that the choice of Vg do 
not alter the allowed range of parameter space. Observation of Fig. yields that, apart from SM 
Higgs resonance region {m ~ mi 12) there exists another allowed range of m — cisi parameter space 
in the vicinity of non-SM scalar resonance (m ~ m 2 / 2 ). Note that variation of m with usi depicted 
in Fig. depends only on the masses of scalar bosons and does not suffer any significant change 
due to change in Vg. The non-SM Higgs signal strength R 2 (calculated using Eq. 16) for the valid 
fn — fJsi parameter space shown in Figs, [^-b is very small and R 2 < 0.2. 

In this work, we assumed two values for VEV Vg (246 GeV and 500 GeV) for the hidden sector 
Higgs doublet <hHS- Erom Eq.[^ we observe that the mixing between the scalars hi and /i 2 depends 
on the VEV of <I>hs and H. Hence, the choice of Vg may change the range of available model 
parameter space. In Eigs. [^-b, we plot the variation of Higgs mixing angle a between hi and /12 
with A 3 for 7712=100 GeV and 110 GeV with mi=125 GeV (mass of SM-like Higgs). Needless to 
mention the region of a — A 3 space shown in Eigs. [^-b are consistent with the bounds form vacuum 
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(a) (b) 

Figure 3: The valid model parameter space in A 3 —a (in deg) plane obtained for the case of m 2 = 100 
GeV (lef panel) and m 2 = 110 Gev (right panel) satisfying the limits from vacuum stability, LHG 
findings, PLANGK DM relic abundance and direct detection limits on crsi from LUX experiment. 

stability, SM Higgs signal strength from LHG, relic abundance of DM from PLANGK and limits 
on DM-nucleon scattering cross-section from LUX direct DM search experiment. Plots in Fig. 
are produced using similar method we have applied previously to obtain viable model parameter 
space for Fig. However a — A 3 plane in Fig. is further constrained by imposing LUX DM 
direct detection bound. The plots in Fig. are for the case when mi = 125 GeV and m 2 = 100 
GeV while plots in Fig. represent the allowed a — A 3 parameter space when m 2 = 110 GeV for 
the fixed value of mi = 125 GeV. The green and blue regions in Fig. and Fig. |3> correspond to 
two different values for VEV of dark Higgs doublet, Vg = 246 GeV and Vg = 500 GeV respectively. 
From Fig. (m 2 = 100 GeV case) one observes that for both the considered values of VEV Vg, 
the mixing parameter A 3 remains small and is confined within the region IA 3 I < 0.01. Eor the case 
when Vg =246 GeV (the red region of Eig. [^), the limit of mixing angle a ranges between —0.1 
to 0.1. However these range (of mixing angle) varies within the limit |a| < 0.2 when Us=500 GeV 
is chosen (green region shown in Eig. [^). Study of the A 3 — a plots in Eig. (plotted for m 2 = 
110 GeV) reveals that for both the values of Vg considered in Eig. the mixing parameter is small 
(IA3I < 0.01). The mixing angle a is bounded in the range |a;| < 0.15 and |a;| < 0.30 for Vg = 246 
GeV and 500 GeV respectively. 
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BP 

Vs 

in GeV 

m 2 

in GeV 

m 

in GeV 

5 

72 

(^si 

in cm^ 

cm^/s 

BPl 

246.0 

100.242 

50.0 

-4.86e-03 

0.60e-06 

2.89e-46 

1.98e-26 

BP2 

500.0 

110.321 

55.0 

-5.85e-03 

0.43e-06 

1.13e-46 

1.90e-26 


Table 1: Benchmark points obtained from the constrained model parameter space in agreement 
with the bounds from vacuum stability, SM Higgs signal strength from LHC, DM relic density from 
PLANCK and LUX DM serach bounds on DM-nucleon scattering cross-section. 

4.2 Calculation of gamma ray signals from galactic centre and dwarf galaxies 



Figure 4: Comparison of the GC 7 -ray flux data from | 22 j with those calculated for benchmark 
points in Table [TJ 

In this Section, we calculate the 7 -ray flux from the galactic centre and dwarf galaxies for 
the fermionic dark matter in framework of the present model and compare our results with the 
experimental observations. For these calculations we consider two benchmark points (BPs) from the 
restricted parameter space that satisfy both theoretical and experimental bounds (mainly vacuum 
stability, LHC constraints on SM Higgs signal, PLANCK results for relic abundance and direct 
detection limit on m — fisi from LUX) for two choices of /12 mass, mainly, 7712 = 100 GeV and 200 
GeV. In Table we tabulate the chosen BPs along with model parameters. There are two chosen 
sets of benchmark points in Table and we denote them as BPl and BP2. The GG gamma ray 
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flux is calculated using Eqs. 31 33 for the BPs tabulated in Table The annihilation cross-section 
{^^)bb dark matter particle is calculated using Breit-Wigner enhancment technique using 

Eqs. |37||^ discussed in Sec. The gamma ray spectrum ^ in Eq. [^is obtained from Ref. m 
for annihilation of DM into any specific channel. The gamma ray spectra for BPl and BP2 are 
then calculated for the specified region of interest adopted from Ref. [22] (|/| < 20*^, 2^ < \h\ < 20°) 
using NEW halo profile (with 7 = 1.2, pq = 0.4 GeV cm“^). In Eig.|^ we show the calculated GC 
gamma ray flux (in GeV cm“^ sr“^) for our proposed DM candidate with BPl and BP2. We also 
show in Pig.j^the CCW data for comparison. Green and blue lines in Fig. [^represent the calculated 
7 -ray spectra for BPl and BP2 respectively. Both the benchmarks points are in agreement with 
the findings from GC gamma ray study presented in CCW [22j. From Fig. it can be observed 
that flux calculated using the set BPl (m=50 GeV) is in better agreement with the hndings from 
CCW analysis. 

We now further investigate how well the DM candidate in our model can explain the observed 
extragalactic 7 -ray signatures from various dwarf galaxies. From their six years observations on 15 
dwarf galaxies, the Fermi-LAT experiment did not obtain any significant excess of 7 -rays. Fermi- 
LAT collaboration m however in a recent work provides combined bound on DM mass and ther¬ 
mally averaged DM annihilation cross-section into SM particles for these 15 dSphs. A similar bound 
in m — (crv)j plane is also presented recently in an another work (T^j for eight new dSphs jointly 
by Fermi-LAT and DES collaboration. In this work we calculate thermally averaged annihilation 
cross-section of DM annihilating into SM sector in our model and compare them with experimental 
results given by [ITl[l2|. In Fig. we plot the bounds on DM annihilation cross-section (for 

the annihilation channel DMDM—)• bb) with dark matter mass m obtained from galactic |25| and 
extragalactic miiiii 7 -ray search experiments. We calculate the variations of the same plotted in 
Fig.j^for the benchmark points BPl (for m 2 = 100 GeV) and BP2 (for m 2 = 110 GeV) considered 
in our model. 

Black contours shown in Fig.j^are the la, 2a and 3a contours given by the CCW [22| analysis of 
GC gamma ray excess observations. The blue line in Fig. [^describes the bounds in m— {aY)p^ plane 
given by the analysis of gamma rays from previously discovered 15 dSphs and they are adopted 
from m- Also shown in Fig. the yellow band which is the 95% conhdence limit (C.L.) region 
adopted from the analysis in Ref. |llj for DM annihilation into bh. The combined bounds on (cv)^j 
for different DM mass m from a recent study of the newly discovered 8 DES dwarf galaxies [T2| are 
given by the pink coloured line in Fig. The green horizontal line in Fig. shows the annihilation 
cross-section for thermal dark matter that may yield the right DM relic abundance obtained from 
the PLANCK experiment. 

From Fig. [^one readily observes that the calculated values of (o'v)^^ for the benchmark points 
BPl and BP2 in our model broadly agrees with the la, 2a and 3a allowed regions in m— (cv)plane 
obtained from the experimental results. This can also be noted from Fig. that these benchmark 
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Figure 5: The allowed range of m — (uv) space (for annihilation into bb) along with the bounds on 
(crv) (into bh channel only) obtained from GC 7 ray search results CCW [ 22 ] and dwarf galaxies 
[111 [T 2 | compared with the same obtained from benchmark points in Table 


points are consistent the combined limit from DES dwarf satellite data and falls within the 95% C.L. 
limit predicted by Fermi-LAT for 15 dSphs. Also the calculated values of for the benchmark 

points considered in our work lie below the upper bound on thermal DM annihilation cross-section. 
Hence, DM fermion in the present model can account for the galactic centre excess in 7 -ray and is 
also consistent with the bounds on gamma ray flux from Milky-Way dwarf satellite galaxies. 

We now calculate the gamma ray flux for 8 new dwarf satellite galaxies discovered by the DES 
experiment for the hidden sector fermionic dark matter candidate proposed in this work. These 
calculations are performed with each of the benchmark parameter sets BPl and BP2 given in 
Table The Gamma ray flux for each of these 8 dSphs in the work |12j is computed using Eq. 34 


and the values of the J factors (Eq. 35) for each of the eight dSphs adopted from Ref. |T2|. In 
Ref. [T 2 | these J factors are estimated by integrating the dark matter density (adopting NFW halo 


profile for DM density distribution) along the line of sight over a solid angle AD = 2.4 x 10 
sr“^. As previously mentioned the gamma ray spectrum ^ is also obtained from Ref. m for 
this calculation. The calculated flux for each of the eight dSphs are shown in eight plots (a-h) of 
Eig.§ Also shown in each of the eight plots of Fig. [^ the respective upper bounds of the flux 
given by the experimental observations of gamma rays from each of the eight dSphs. These are 
shown as red coloured points while the computed flux in this work for the respective dSphs are 
given by continuous lines in Fig. The green and blue continuous lines in each of the plots (a-h) 
of Fig. [^correspond to the calculated flux using the benchmark points BPl and BP2 respectively. 
It is clear from Fig. that the fluxes calculated, assuming the annihilation of the DM candidate 
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(b) 


(c) 





(d) (e) (f) 




(g) (h) 

Figure 6 : Comparison of the observed upper bound on 7 -ray flux for 8 DES dSphs with the 
calculated 7 -ray flux from BPl and BP2 tabulated in Table 
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in our proposed model, for all the eight dSphs do not exceed the upper limit of 7 flux set by the 
experimental observations of DES collaboration. 



Figure 7: Benchmark points BPl and BP2 compared with the allowed region of model parameter 
space shown in m — Ji 9 ((tv )_26 plane obtained from [TO] 

Besides the 15 dwarf galaxies investigated earlier and the eight other recently explored dwarf 
galaxies, one more dwarf galaxy namely Reticulum 2 (Ret2) has been probed very recently. Geringer- 
Sameth et. al. uni, after an analysis of observed gamma rays from Ret2 dwarf galaxy reported 
an excess of gamma ray emission from Ret2. From their analysis of Ret2 data Geringer-Sameth 
et. al. provide different C.L. allowed contours in m — Jig{av )-26 plane where m is the mass of 
the dark matter and Ji 9 (crv )_26 is the product of the J factor in the units of 10^® GeV^ cm“® and 
thermal averaged product {av) of annihilation cross-section and relative velocity in the units of 
10“^®cm^s“^ for various final state SM channels. As mentioned earlier in this work DM candidate 
primarily annihilates into bb, only the contours for the DM pair annihilation into bb channel are 
adopted. For the present dark matter model with the constrained parameter space discussed earlier 
we compute the quantity Jig{av)- 2 e for different dark matter mass m annihilating into bb channel. 
However the value of the J factor for Ret2 has been adopted from [TO]. In their work Geringer- 
Sameth et. al. |in) estimated the J values by performing line of sight integral over a circular region 
with angular radius 0.5*^ surrounding the dwarf and over a solid angle AD = 2.4 x 10“^ sr“^. All 
these calculations are performed for two values of non-SM scalar mass accounted in the present 
model namely m 2 = 100 GeV and m 2 = 110 GeV. The results are presented for the two benchmark 
points BPl and BP2 corresponding to the calculations with m 2 = 100 GeV and m 2 = 110 GeV 
are shown in red and skyblue points in Fig. In Fig. the contours from the experimental data 
analysis by Geringer-Sameth et. al. are given for comparison. In Fig. the contours for 68 %, 


24 







95% and 99.7% C.L. are shown in black coloured lines in increasing order of area enclosed by each 
contour. The valid regions of m — Ji 9 ((Tv )_26 plane in our model (calculated for DM annihilating 
into bb pair) are presented by green coloured patches in both the plots of Fig. From Fig. it 
can be easily observed that Ji 9 ((Tv )_26 in the present model calculated for DM annihilating into 
bh channel (for benchmark points with m2=100 GeV and 110 GeV) is within the 3a C.L. limit. 
Hence fermionic DM candidate in the present framework can also explain the observed excess in 
7 -ray from Ret 2. 

5 Discussions and Conclusions 

In this work, we have proposed the existence of a hidden sector which obeys a local SU(2)h and a 
global U(1)h gauge symmetries. In order to introduce fermions which are charged under this SU(2)h 
gauge group one should have at least two fermion doublets in order to avoid “Witten anomaly”. 
The particle and the antiparticle of these dark fermions are different as they possess equal and 
opposite U(1)h charges. Similar to the usual Higgs doublet in the visible sector, this hidden sector 
also has an SU(2)h scalar doublet which however does not have any U(1)h charge. The SU(2)h 
gauge symmetry breaks spontaneously when the neutral component of the scalar doublet gets 
a VEV and thereby generates masses of all the dark gauge bosons (%(j) and dark fermions (/j). 
Since the dark sector fermions interact among themselves through the dark gauge bosons, therefore 
all the heavier fermions as well as the dark gauge bosons can decay into the lightest fermion and 
hence the lightest fermion in this dark sector can be treated as a particle for the viable dark matter 
candidate. In fact in this model this lightest fermion is the only dark matter candidate. The dark 
fermions and dark gauge bosons do not mix with the SM fermions and gauge bosons due to the non 
abelian nature of two SU(2) groups. However, the dark sector scalar field can interact with the SM 
Higgs like scalar in the visible sector and only through this interaction two sectors are mutually 
connected. 

We therefore test the viability of the present model by using theoretical and experimental 
constraints on the relevant model parameters, such as vacuum stability conditions, bounds on relic 
abundance of DM from PLANCK experiment, direct detection limits on DM-nucleon scattering 
cross-section from LUX experiment. LHC bounds on signal strength and invisible decay width 
of the SM Higgs, are also used to constrain the parameter space. From such analyses we find 
that only a small region of the parameter space near the scalar resonances (when m ~ ^ and 
^), is consistent with the current experimental bounds. Study of the model parameters, thus 
constrained, shows that the mixing between the two scalars {hi, / 12 ) of the model is very small 
(mixing angle a < 0.3 deg) and depends on the VEV of the dark scalar doublet (^h). With the 
allowed regions of parameter space, thus obtained, for the present DM candidate (dark fermion /i) 
we compute the gamma ray flux from the GC region. While calculating the gamma ray flux from 
GC we have used Breit-Wigner enhancement mechanism for the computation of DM annihilation 
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cross-section into the bb final state (/i/i —)• bb). These computational results are then compared 
with the experimental analyses of the Fermi-LAT GC gamma ray flux data considering the dark 
matter at the GC primarily pair annihilates into bb channel. Our proposed DM candidate can 
indeed explain the results from these experimental analyses. 

In search of indirect evidence of dark matter from astrophysical sources, the gamma rays from 
various dwarf satellite galaxies are also explored for possible signature of excess gamma rays from 
these sites. To this end 15 such dwarf galaxies have earlier been investigated and more recently 
the gamma ray observation is also reported from eight more newly discovered dSphs. From the 
analyses of these observational results different C.L. bounds have been given in the parameter 
space of (o'v)^j; — m plane. We compare our computational results with these experimental bounds 
and found that the 7 -rays that the DM candidate in our model produce on pair annihilation can 
simultaneously satisfy the observational results from GC and dwarf galaxies. We also demonstrate 
that the calculated fluxes in our model for each of the recently discovered eight dwarf galaxies lie 
below the corresponding upper limits of the fluxes obtained from the observational results of these 
dwarf galaxies. We further demonstrate that our calculations are also in good agreement with the 
analysis of Ret2 dwarf galaxy observations. 

Our work clearly demonstrate that the dark matter candidate proposed in this work is a viable 
one to explain the 7 -rays from both the GC region and dwarf galaxies simultaneously. However 
the dark matter can also pair annihilate into fermion-antifermion pairs and there are experiments 
such as AMS-02 that look for the excess of e^/e“ or pp in cosmos. In a recent work, AMS- 
02 collaboration have reported their first measurement of pjp flux [119] . A model independent 
analysis of this AMS-02 p/p data is performed by Jin et.al. [120j . In this work |120| . the upper 
limits in (crv) value for DM annihilation into SM particles (quarks and gauge bosons) for different 
considered DM halo profiles (NFW, Isothermal, Moore) are obtained. The analysis presented in 
the work |120] also considered four different propagation models namely conventional, MED, MIN 
and MAN[^ We have also checked that the DM in our model satisfies upper bound on (o'v)^^ given 
in Ref. |120j when NFW profile is considered. This is found to be true for both the cases of dark 
sector scalar mass m 2 = 100 GeV and m 2 = 110 GeV. Hence, fermionic dark matter explored in 
the present model can serve as a potential candidate for dark matter. Upcoming results from LHC 
as also DM direct and indirect search experiments may provide stringent limits on the available 
model parameter space. 
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Appendix A 

®For further studies see mi and references therein. 
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Coupling between the scalars hi and /i 2 are given as follows 
Am = Xivcl - X2VsSa + ^Xsivcasl - VsSacl) , 

A 222 = + A2 Us4 + ^hivSacl + VsCaS^) , 

A 211 = 3{Xivclsa - X 2 Vsslca) + ^A3(r;s(cj^ - 2s^Ca) + v{s^ - 2clsa )), 

A 122 = 3(Air;S^Ca - X 2 Vsc\sa) + ^A3(r;s(-4 + 2c\sa) + v{c% - 2s\ca)) . 
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